Thermodynamics of boson and fermion systems with fractal distribution functions 
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Starting with the fractal inspired distribution functions for Maxwell-Boltzmann, Bose-Einstein 
and Fermi systems, as reported by F. Biiyiikkilig and D. Demirhan, we obtain the corresponding 
probabihty distributions and study their thermodynamic behavior. We compare our results with 
, those corresponding to ideal gases {q — 1), and Bose-Einstein and Fermi systems with quantum group 

' symmetry. In particular, we show that the hamiltonian that gives the Bose-Einstein generalized 

^\ I distribution function can be interpreted as a g-deformation of the ideal gas hamiltonian. 

U , 

Ph, PACS number(s):05.30.-d 
< 

^ . I. INTRODUCTION 

> 

' As is well known from the theory of fractals, given the statistical weight Q{q, r) of a system and the resolution r, the 

» . . . ... 

, fractal dimension is defined as the exponent d = Dq which will render the product limr^o ^{q, f)T finite. Since the 

■ entropy S{q,r) is proportional to lnQ{q,r), the relation between the entropy S{q,r) and the generalized dimensions 

o^ : 

Q\ Dq IS given by 



2: ^- lim^'^) 

r^o Inr 



X 



(— I ' where the q parameter stands for the moment order. The entropy function will depend on a set of probabilities pi that 
O 

^ J a random variable will be found into the ith bin of size r. Based on these definitions and by use of the Boltzmann's 



_ff-theorem, the generalized entropy and distribution functions for classical and quantum gases were calculated in Ref 
101 . The average number of particles with energy e were shown to be given by 

where a = for the classical case and the values a = — 1,+1 correspond to Bose-Einstein and Fermi-Dirac cases 
respectively. For 9=1, Eq.(H) becomes the standard textbook result for classical and quantum gases. Therefore, 
for arbitrary q and within the context of statistical mechanics, the corresponding formalism can be understood as 
a generalization of Boltzmann-Gibbs statistics. Two examples of generalizations of the concept of entropy based on 
ideas from the theory of fractals are the functions introduced by Renyi |^ and Tsallis , which are given respectively 
as 
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where pft. is the probabihty for the ensemble to be in the state R. It is clear that both entropy functions become the 
Shannon entropy function S = —k "^j^Pr Inpij as q ^ 1. A discussion on the irreversible character of these entropies 
can be found in Refs. Q|. In particular, the entropy function Sq has all the properties of the Shannon entropy except 
that of additivity. In fact, given two independent systems S and E' the function Sq satisfies 

which becomes additive for q — 1. Therefore, Tsallis' formulation provides a framework to deal with the nonextensive 
properties of certain physical systems. Numerous applications of this formalism 1^, and similar formulations appeared 
in the literature . The equilibrium probability distribution pi^ in the grand canonical ensemble is written as 

[l + fi{q-l){ER-^xN)f"^'~'^^ 
PR = 7^ > (6) 



with the partition function 



Zq = J2[l + I3{q l){Eii ^,N)f^'-'^'^ , (7) 

R. 

and the total energy Er — njSj. Thermodynamic stability is achieved by defining a q- averaging 

{E)^Y.Pr^R^ (8) 

R 

leading to expressions of the type 

d Z]-'i - 1 

The fact that the partition function in Eq. does not factorize in independent modes is a consequence of the nonex- 
tensive characteristics, like those involving long range interactions, of a physical system. Due to the mathematical 
complexity of dealing with equations like (^, calculations in this direction have been performed |^ for values of (7 sa 1. 

The distribution functions in Eq. and nonextensive thermodynamics are not unrelated. In fact, these fractal 
distributions have also been obtained |^ by considering the case of a dilute gas and approximating the partition 
function in Eq. (Q) with a factorized partition function. This approximation to the nonextensive thermodynamics of 
Ref. which ignore the correlations between the particles, has been shown [|o| to be a good one outside a certain 
temperature interval. This time interval shifts to higher values of T when either the number of particles or the 
number of energy levels are increased. Some applications of the dilute gas approximation as an approximate approach 
to nonextensive thermodynamics can be found in Refs. 
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It has been also pointed out ||T2[ that nonextensive thermodynamics could be understood in terms of ^-deformations. 
The fact behind this possible connection between two apparently unrelated subjects resides basically in the observation 
that a (/-number [x] = {q^ — l)/((7 — 1) has the pseudoadditivity property of the type displayed in Eq. and thus 
the entropy Sq can be defined in terms of the g-derivative — x^^{l — — q) acting on the probability 

distribution. The theory of q-analysis was formulated at the beginning of this century and it serves as an analytical 
framework to study basic hypergeometric series |Q . The study of the so-called g-series dates back to the times of Euler, 
and it has played an important role in the theory of partitions [ p^ and more recently in quantum group theory [ p^ . 
Although nonextensive thermodynamics does not embody any quantum group structure, it could be somehow related 
to the theory of quantum groups. The link to quantum groups is through (/-analysis. As is well known, noncommutative 
calculus arises as a set of algebraic relations which are covariant under quantum group transformations [ p^ . It has 



been shown |17| that differential operators in noncommutative calculus has a representation in classical space in terms 
of (/-derivatives times scaling operators. Therefore, although the link between nonextensive thermodynamics and 
quantum groups seems to be rather marginal, these two distinct subjects share a common language which is that 
of (/-analysis. In addition, it has been shown that generators of multifractal sets which are self-similar under 
(/-scalings can be identified with the (/-derivatives Dx- 

Then, according to the work of several authors, the fractal distribution functions in Eq. are those corresponding 
to a dilute gas approximation of nonextensive thermodynamics, and the formalism of Tsallis thermodynamics seems 
to be related to the theory of quantum groups through the use of (/-analysis. Therefore, it is natural to investigate 
whether the thermodynamic functions arising from the fractal distribution functions in Eq.(|^) are somehow similar 
to the thermodynamic behavior of systems with quantum group symmetry. The meaning of the parameter q should 
be understood within the appropriate context. In the theory of fractals q stands for the moment order that defines 
the fractal measure |l9| ], in the theory of quantum groups q parametrizes a deformation from classical Lie group 
symmetry, and in Tsallis thermodynamics q parametrizes the departure from extensive thermodynamics. In this 
work, the most appropriate interpretation of (/ is as a parameter that measures the departure of the corresponding 
statistical mechanics from the Boltzmann-Gibbs formulation. 

Our motivation in studying this system is twofold. First, as discussed above, it relates to the interesting fact 
that the fractal distributions functions approximate well, within certain limits, nonextensive behavior. Second, it 
is of theoretical interest to study a thermodynamic system whose probability density generalizes Boltzmann-Gibbs 
statistics. Our calculations will also indicate whether the corresponding behavior from these generalized statistics 
shares some properties with the thermodynamics of quantum gases with quantum group symmetry. 

As our starting point we consider the two following basic requirements. In order that the applicability of the 
thermodynamic limit makes sense we consider a factorized probability density. The formulation should also be 
consistent with the classical limit. 
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This paper is organized as follows. In Section 2 we obtain the probability density and density operators for Maxwell- 
Boltzmann and quantum statistics respectively that lead to the particle distributions functions in Eq.(|^). For the 
case of Maxwell-Boltzmann statistics, we find the correct way of calculating the thermodynamic functions such that 
the expected classical behavior occurs. For quantum statistics, with use of the standard operator formalism we derive 
the boson and fermion Hamiltonians that leads to the occupation number functions in Eq.(^, and based on this we 
study the thermodynamics of these systems in the thermodynamic limit. In particular, for the Bose-Einstein case we 
calculate the dependence of the internal energy, heat capacity and entropy on the parameter q and the temperature. 
For quantum statistics we compare our results with the thermodynamics of quantum group invariant systems as 



reported in Refs. |20 2l|]. In Section 3 we summarize our results. 



II. THERMODYNAMICS OF BOSON AND FERMION SYSTEMS 
A. Maxwell-Boltzmann case 

The particle distribution for the Maxwell-Boltzmann case is given by Eq. for a = 0. 

{ni)^pr\ (10) 

where 

= (1 + (q - l)/3(e; - . (11) 

We define the probability density p according to the equation 

with the corresponding partition function 



l=Qni=0 

-1 



He"'" . (13) 



1=0 

A simple check shows that from the definition of the average total number of particles {N) 

En,^0(l/»i!)^/P'"' 

E„,=o(l/^^!)p""' ' 



-E.-^^> (-) 

1=0 ^ 

the probability density in Eq. ( p^ leads to the particle distribution in Eq. ( |l0[) . In the thermodynamic limit we rewrite 
Eq. ( p^ ) in terms of the integral 

4:TtV / 2m ^ 



U(9-l)j Jo il + x^-{q^lWf^'-'y ^^^^ 
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For further purposes we wish to consider the integral 

Is = 



(16) 



(l + a;2-((7-l)/3/i)i/(«^^'' 
The requirement of finiteness for Is restricts the values of q to the interval q G [1, (4 + 5)/(2 + (5)]. Defining a new 
variable w = 1 + .x^ — (g — and expanding the integrand in powers of ((g — — 1) /w we find that Is is given 
by the expression 



1 



where Ss is a series independent of Pfi given by 

1 (-1)™ fS-l\ fS + 1 



l)/2 



Ss 



E 



-l/(g-l) + (<5 + l)/2 .... 
Therefore, for Eg. (p^) we write 

-2TrV f 2m 



1 



-l/(g-l)-m+(5 + l)/2 



{N) = 



3/2 



\Piq-l)J (1 - (g - l);3Ai)l/('-l)-(3/2) 

such that the temperature dependence of the fugacity z is given by 



(17) 



(18) 



(19) 



In z 



9-1 



2TrV f 2m 



3/2 \ 2(g-l)/(5-3g)- 
^2 



(20) 



restricting the fugacity to the interval < z < e^/^'"-'^^ An inspection of Eq.(pT[) shows that the correct way of 
defining the average energy is through the equation 

AirV p p2 



/l3 



(21) 



such that the average energy per particle becomes 

([/) _ kT J^x\l + x^)-^^(^-^Ux 



(N) Q-l/g°°x2(l + x2)-i/(9-i)rfa;' 
Convergence of these integrals restricts the values of q to the interval 1 < g < 1.5 With use of the simple identity 



(22) 



dx = 0, 



/o dx(^i + ^2fii-^) 
it is easily checked that Eq.(|2|) leads to the required classical result 



([/) - -{N)kT. 



(23) 



(24) 



Therefore Eq. (pi|) indicate how we should calculate the average energy for boson and fermion systems such that the 
classical behavior be recovered in the limit of high temperatures. 

In order to calculate some thermodynamic properties of boson and fermion systems with particle occupation numbers 
given by Eq.(||), we consider the general density operator 



'^=|n(i+/'(9-i)^0'^^' (25) 

with the partition function Z given by the expression 

^ = n E (i + - 1)^^) '^^'"'^ ' (26) 

j=0 nj=0 

such that Trp = 1. We wish to find the form of the function K — Kj which would lead to the average occupation 
numbers given in Eq (|^) and the standard thermodynamics expressions as q — s- 1. It is natural to assume that 
the function will depend on the energy e^, the chemical potential /i and the number operator hj. Following our 
procedure for the classical case, we will calculate the occupation number and the average energy from the corresponding 
definitions (iV) = TrpN and (t/) = TrpiR. 

B. Boson case 

With use of the standard boson algebra 

afc'a|. - aJ,Ofc' = 5fe,fc', (27) 

and assuming that aJ,/5 = pkpal, it is simple to check that the expectation value of the product of a creation and 
annihilation operator can be written as 

{alak') = -SkM' + Pk{a\ak'), (28) 

leading to 

(4«^'> = (29) 
Pfe - 1 

where the function pk is given by Eq.(^). Based on the observation that for an arbitrary function (7, independent of 
the operators a and a\ and the number operator fij = ajoj we have the identity g^^a^^ — ga^g"^, we find that the 
function Kj is given by the expression 

where we need to restrict q to those values such that l/{q — 1) is an integer. Therefore, the parameter q = {j + l)/j, 
where j is a positive integer. For those values of q such that l/(g — 1) is an odd number the function n(e) > for 
all values of the energy and p < 0. For l/{q~ 1) even, although n{0) > for /3/i ^ [0,2/(5 — 1)], the derivative 
d{n) jde < for /i < only. Therefore, as in the standard Bose-Einstein case the chemical potential is negative for 
all temperatures T > 0. 

The density operator becomes 
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with the partition function Z given by 

Zbe = n T^-i- (32) 



Equation ( pO| ) can be rewritten in a more suggestive form by defining the variable Qu — p\ ^ such that the full 
operator K reads 



k^Y.^e,~y.)\^, (33) 



which clearly contains a sum of products of the factor (sj — /i) times a Qj-number. As a matter of fact, we can rewrite 
Eq. ( |3^ ) in a more standard form as follows 

i: = ^(e,-M)0,0,, (34) 
where the set of operators and its adjoint are defined in terms of boson operators according to 

!>, = a] (35) 

*,=«;-^. (36) 

The operators 0j and 0^ satisfy an algebra reminiscent of the so-called g-boson algebra, which is 

(^kf^j - Qk'''4>j4'k = Sj^k- (37) 
For the average internal energy (U) ~ Trp'^Hs, with Hb ~ ^k4'k4'k, we simply get 

k Pk-Pk 

In order to evaluate Eqs. (^) and ( |38| ) in a simple way we approximate these equations by considering the thermody- 
namic limit. Therefore we write 

{N)^{N{0)) + ^Gy2iz,q), (39) 

(C^) -^G5/2(z,g), (40) 

where we have separated, as usual, the ground state occupation number, — h? /2'KmkT , and the functions Gu are 
written 

4 /■°° 

G:i/2{z.q)^ —j=-, 7^372 / T .i/(a-i^ (^1) 

8 [°° 

^^/^(^''^)=3V^(,-l)a/^yo [l + f7r/(«-)-(l + 0)'"' ^''^ 
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where fl 



{q — l)lnz. The functions 0^/2(2, q) and 0^/2(2, q) become for q = I the well known Bose-Einstein 



functions = J2n=i ^"Z"-" for v = 3/2,5/2 respectively. Figure 1 shows the graphs of the functions 63/2(1,'?) and 
6*5/2(1, 9) in comparison to the functions 32/3(1, ?) and 55/2(1, q) for the quantum group case in Ref. pl[ |. As expected 
these functions converge at 5 = 1 to the values C(3/2) = 2.612 and ({5/2) = 1.341 respectively. These graphs indicate 
that the thermodynamic behavior for the fractal case will be different than the case of quantum group bosons. In 
particular, since the critical temperature Tc is proportional to 1/63/2(1,'?)^^'^, we see from Figure 1 that the critical 
temperatures for a Bose gas and for a quantum group boson gas T^'~^ compare to Tc as follows 



Defining the function (u) = {2T:m/h^){V^ / {N)^y/^ (U) in Eq.(0 

, > 3 65/2(2,9) /T 



^Gl%il,q)\Tc 



we write 

5/2 



(43) 



(44) 



Fig. 2 shows the function {u){T) for q = 1.1 and q = 1.2, for an ideal gas and quantum group SUi,i{2). From the 
figure we see that U'^^ > U^^ > U and solely the quantum group case exhibits a discontinuous derivative at the 
critical temperature. The heat capacity Cy = {d(U)dT)y is given by 

"1565/2(^,9) 963/2(2,9)^^^ 



CviT) = k{N) 



463/2(2,9) 463/2(1,9) dGsMz.,) 

oz 



3/2 



which reduces for T < Tc, /i = and dfi/dT = to 



Cv{T<Tc)^^{N)k(^ 



3/2 



65/2(1,9) 



(45) 



(46) 



63/2(1,9)' 

Since the function dG^/2/dz is not finite in the z ^ 1 limit, there is no gap in the heat capacity at T — Tc- Figure 3 
shows a graph which compares the heat capacity for 9 = 1.1, as given by Eq.(^), with the textbook case of an ideal 
Bose-Einstein gas and the quantum group case as reported in Ref. [ pl| . At high temperatures, the heat capacity Cy 
in Eq.(|4|) tends to the classical result 6y = {3/2)k{N). The entropy function S = ((1//3) InZ + {U)) - n{N)){l/T) 
is given by the expression 



S ■ 



k(N) 



63/2(1,9) 



365/2(2,9) 



(9 



.^ln(l-(l + .)Va-.)),.-ln.] (I) 



3/2 



(47) 



Figure 4 shows a graph of the entropy per particle for 9 = 1.1 in comparison to those corresponding to the ideal Bose 
gas and the SUi,i{2) cases. For the quantum group SUi,i{2) the entropy is given by the expression 



5 35/2(^,9) / T 

2 ^ ^3/2(1,9) Ue^^ 



3/2 



In; 



(48) 



From these graphs we see that the behavior near the corresponding critical temperature is very different than the 
quantum group case in which the heat capacity exhibits a gap that increases with the value of the parameter 9 and 
the entropy has a discontinuous derivative. 



C. Fermion case 



For the fermion case it is simple to check that 

n (i + p{q - m) bi = (1 + p{q - i){ek - i^))bi n (i + P{q - mj 



(49) 

where the operator Fj = [sj — ii)b%j. Therefore the occupation nmnber {b%j) = Trpb%j and the average internal 



energy ([/) — Trp'^Hp become 



(l + ;J(«-l)(e, -,,))■/(«-■) 



(l + ;J(<)-l)(e, 

J 

We consider, as in the previous section, the thermodynamic limit of Eqs.(pO|) and ( ]5l| ) 

W / 1 



(iV) = 



P3/2iz,q), 



W I 1 1 



1)5/2 ^3^- 



where the functions Fiy{z, q) are defined by 



F- 



3/2 



F5/2 — 



1 + [1 + f7] 

^4 



1/(9-1) 



-dx. 



lo 1 + [1 + r!]«/(«-i) 

Performing a similar calculation already done in the previous section, we find for the heat capacity Cy 



15 



dz 
_(f,9) 
dz 



kT 



3/2 



(50) 
(51) 

(52) 
(53) 

(54) 
(55) 

(56) 



- (^3T)i75 y^5/2(^,g) - -F,/,{z,q)^^ 

\ dz 

where p\l^ = i3{N) /4TrV){h^ /2mf/^ is the Fermi energy. For the quantum group SUq(2), the heat capacity is given 
by the equation 

15 9 9h/2iz,q) ■ 

where the functions /3/2 and are given as follows 



C^^ = k{N) 



dz 
l/2{z,q) 
dz 



kT 



3/2 



(57) 



fs./2{z,q) = 



h/2{z.q) 



00 X' 



dx. 



1 

x^ (2 + (q-2 + i)e-9-'-') 
1 



dx. 



(58) 



(59) 



2 1 —2 2 

with / = +2 + er'^ ^ z. Figure 5 shows a graph of the heat capacity for q = 1.1 and q = 1.2, the ideal Fermi 

gas and the quantum group SUi,i{2) fermion case as a function of the temperature. In general, the graphs show that 
the heat capacities follows the relation 



Cv < < C^^ 



(60) 
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III. CONCLUSIONS 



In this paper we have studied the thermodynamics of the fractal distribution functions for MaxweU-Boltzmann, Bose- 
Einstein and Fermi systems in the thermodynamic hmit. As pointed out in the Introduction, the parameter q measures 
the departure from the Bohzmann-Gibbs statistics. As our starting point, we require that the formahsm be based 
on a factorized probabihty density and be consistent with the classical limit. Based on this, with use of the operator 
formalism we obtained the matrix density operators that lead to the quantum generalized particle distributions 
n{e). For the Bose-Einstein case, we obtained the interesting result that the hamiltonian can be interpreted as a 
q-deformation of the boson ideal gas hamiltonian. We should remark that for the boson case our partition function 
differs from the factorized partition function considered in Ref . |^ . A calculation of the heat capacity that results from 
the boson generalized distribution function shows that it is continuous at the critical temperature. This contrasts 
with the quantum group case wherein the heat capacity exhibits a discontinuity at that temperature that increases 
with the value of the deformation parameter. For the fermion case, the heat capacity function for the Fermi ideal gas 
lies between the heat capacity for the fractal and quantum group cases. In general, our results point out that the 
thermodynamic behavior given by both fractal distribution functions are very different to those found for Bosc-Einstcin 
and Fermi systems with quantum group symmetry. Since, as claimed for several authors, the fractal Bosc-Einstcin 
and Fermi distribution functions approximate well, under some conditions, the nonextensive behavior of some physical 
systems, the relation that the latter may have with quantum group symmetric systems seems to be rather limited to 
the context of g-analysis. 
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FIG. 1. The functions 03/2(1,(7) and 05/2(1,9) as defined in the text, and the functions (73/2(1,'?) and .95/2(1,'?) 
from Rcf. for the quantimi group, SUq{2), case for the interval 1 < 9 < 1.25. 

FIG. 2. The function u, as defined in the text, for q = 1.1 and q = 1.2, for an ideal Bose-Einstein gas, and for 
the SUi,i{2) boson case as functions of T/Tcriucai, where Tcriucai refers to the critical temperatures Tc, T^^ and 
j-QG^ respectively. 

FIG. 3. The heat capacity Cv as a function of T/TcriUcai showing the sharp difference between the phase tran- 
sition for the quantum group case in comparison to the Bose-Einstein and fractal cases. 

FIG. 4. The entropy function S for the fractal 9 — 1.1, Bose-Einstein and 5't7i.i(2) cases. 

FIG. 5. The heat capacity as a function of the temperature for q = 1.1 and q — 1.2, the ideal Fermi gas and 
the quantum group SUi,i{2) fermion cases. 
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